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Abstract. In this paper, we address the issue of the stability of the thermal 
equilibrium of large quantum systems with respect to variations of the thermal contact 
between them. We study the Schrodinger time evolution of a free bosonic field in two 
coupled one-dimensional cavities after a sudden change of the contact between the 
cavities. Though the coupling we consider is thermodynamically small, modifying 
it has a considerable impact on the two-point correlation functions of the system. 
We find that they do not return to equilibrium but essentially oscillate with a period 
proportional to the length of the cavities. We compare this coupled cavities system with 
the perfect gas which is described by similar expressions but behaves very differently. 
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1. Introduction 

The physically relevant degrees of freedom of a large isolated system initially out of 
equilibrium, are expected to relax to their thermal equilibrium values. Obviously, no 
relaxation behavior can be observed for the complete quantum state of the system but 
interesting degrees of freedom can evolve irreversibly in the limit of a large system. This 
issue of thermalisation in isolated quantum systems has been principally investigated 
by considering composite systems consisting of two distinguishable subsystems of very 
different sizes. For such systems, the smaller subsystem relaxes to thermal equilibrium 
if the larger one is assumed to be in an equilibrium mixed state pp. Recently, it has 
been shown that a priori thermal averaging is not essential. The small subsystem can 
thermalise whereas the whole isolated system is in a pure quantum state [2]. 

An isolated system of identical particles cannot be divided into distinguishable 
subsystems and the relevant degrees of freedom are, in this case, the n-particle reduced 
distribution functions such as the particle number density. Relaxation behaviors have 
been obtained for both integrable [3l HJ EE E] and nonintegrable systems [8j [9] . In these 
studies, the system is supposed to be initially in the ground state of a given Hamiltonian 
but its subsequent time evolution is governed by a different Hamiltonian. For example, 
a gas is considered to be initially confined in a subvolume of a larger accessible space. 
The single-particle distributions were always found to relax to asymptotic profiles but 
which are, intriguingly, not all well described by an equilibrium ensemble. 

To better apprehend the thermalisation process in isolated many-body systems, it 
is necessary not to be limited to ground states and to be able to treat also finite initial 
temperatures. As we are concerned with isolated systems, such an initial condition 
must be represented by a microcanonical mixed state, according to the postulate of 
equal a priori probabilities, or, equivalently, by a pure state of macroscopically well- 
defined energy [2} [TOl EH]. In Ref.pJ], the time-dependent particle density of a perfect 
quantum gas, initially in such a state, has been determined. It was found that this 
density relaxes from a thermal profile to a non-thermal one. In view of these results, 
one may wonder how two quantum systems, initially independent of each other but at 
the same temperature, evolve if they are brought into thermal contact. Do their relevant 
degrees of freedom remain at equilibrium once they are coupled ? 

In this paper, we address this issue by considering a free bosonic field in two one- 
dimensional cavities coupled to each other by one of their ends. Free fields are commonly 
used to describe the large number of environmental degrees of freedom which give rise 
to dissipation in a smaller system [U [12]. Here, we are interested in two systems of 
comparable size, neither of them acts as a heat bath for the other one. In section [2j 
we present our model Hamiltonian of two quantum systems in thermal contact. This 
model allows to describe from perfectly coupled to completely uncoupled systems. We 
show that the considered coupling between the two cavities is thermodynamically small 
by evaluating the equilibrium total energy and correlation functions. In section [HI we 
determine how the energy contents of the cavities and the field-field correlations evolve 



Sudden change of the thermal contact between two quantum systems 



3 



in response to a sudden change of the coupling between the two cavities. Our model 
is simple enough to obtain exact results for arbitrary values of the coupling strength. 
In section HI we discuss why, though the expressions are formally similar, relaxation 
behaviors are obtained for a perfect gas but not for the system studied here. Finally, in 
the last section, we summarize our results and mention a few open questions raised by 
our work. 



2. Model 



2.1. Derivation of the Hamiltonian 

To describe two coupled one-dimensional bosonic systems characterised by length L and 
velocity c, we start from the Hamiltonian 

,.2 



H'=- I dx 



2 



L 



C 



n(x) 2 



(1) 



where the fields II and (f) are canonically conjugate to each other. The position-dependent 
refractive index n(x) is equal to n > 1 for — d < x < d and to 1 elsewhere. Diverse 
physical systems can be represented by this Hamiltonian. The fields n and d x <p can be 
interpreted as the electric and magnetic components of the electromagnetic field, or as 
the charge and current distributions of an LC transmission line [H[T3]. We write these 
fields as 

U(x) = P + iJlj2 e ~~ 9A/2 ^ q (x) (a{ - a q ) (2) 

' <J>0 

fc* = ^E^ A/2 ^(4+<0 (3) 

where the operators P and a q satisfy the commutation relations [P, a q ] = [a q ,a q i] = 
and [a., aL] = 5 gq i, and A is a cut-off length. Throughout this paper, we use units in 
which h = ks = 1- The eigenmodes ip q of H' are solutions of the differential equation 
~Q 2{ Pq = d x [d x {p q /n(x) 2 ] with the boundary conditions d x ip q (—L) = d x ip q (L) = 0. The 
odd ip q are proportional to sgn(x) cos[g(|x| — L)] for \x\ > d, where the wavenumber q 
obeys 

tan[g(L - d)}' 1 = ntan(2nqd). (4) 

The even eigenmodes are proportional to cos[g(|x| — L)\ for \x\ > d, with q given by 
tan[g(L — d)] = — n tan(2ngd). 

To obtain the Hamiltonian of two cavities coupled to each other by one of their 
ends, we consider the limits n ^> 1 and d-C L with the length 

D = n 2 d (5) 

kept fixed. In this regime, equation (J3|) simplifies to tan(gL) = 2/Dq which can be 
rewritten as 

D « + 2i e -**L = 1 (6) 

Dq-2t [> 
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This second form will be useful in the following. For the even eigenmodes, the 

wavenumbers become q = pir/L where p is a positive integer. It can then be shown, 
using (j2J) and (j3J), that the Hamiltonian (JTJ turns into 

H = i/left + -^coupling + -^right (7) 



where 



Hieft = \f dx [U{xf + c 2 ^) 2 ] , (8) 
-bright is given by this expression with — L and replaced by and L, respectively, and 

^coupling =^ [0(0+) -0(0-)] 2 . (9) 

In the following, we study the Hamiltonian (|7|) which allows to describe from perfectly 
coupled to completely uncoupled cavities. Its eigenmodes ip q obey d x (p g (0~) = 
d x <P q (0 + ) = [<p q (0 + ) — cp q (0~)] /D which leads to even <p q (x) = L" 1 ' 2 cos{qx) 
independent of the coupling characteristic length D. The odd eigenmodes read 

ip q {x) = A q sgn(x) cos[g(|x| - L)\ (10) 

where A q = [L + (D/2)(l + (Dq/2) 2 )- 1 }- 1 / 2 and q is a positive solution of ©. In the 
large D limit, the two cavities decouple from each other as can be seen from the fact 
that the boundary conditions at the ends 0~ and + become similar to that at the ends 
—L and L. In this limit, the solutions q of the equation ([6]) are the multiples of tt/L 
and the corresponding eigenmodes <*p q simplify to L _1 / 2 sgn(:r) cos(gx) for q > and to 
(2L) _1 / 2 sgn(x) for the solution q — » 0. In the opposite limit, D — 0, the odd eigenmodes 
are given by <p q (x) = L -1 / 2 sin(gx) with the wavenumbers q = (p+l/2)ir/L where peN. 
The transmission at x = is perfect in this case. 

2.2. Thermal contact 

The coupling (Q is thermodynamically small, i.e., its contributions to the 
thermodynamic functions of the complete system described by (0) are negligible in 
the thermodynamic limit of large L and fixed energy density. As an example, we 
consider the average energy (H)t = Tr[exp(— H/T)H]/Tr[exp(— H/T)] at temperature 
T. Diagonalizing the Hamiltonian (CO) with the help of the expansions (j2J) and (J3j) gives 



Wt= 2 + |j2tanh(cg/2T)- 

As there is only one solution to equation ([6]) in a given interval \pK/L, (p + l)ir/ L] where 
p EN, the above sum over q can be approximated by an integral, which leads to 

(H) T = + f LT 2 + E (TD/c, TL/c) . (12) 
7iA z 3c 



The first two extensive terms stem from the integral approximation, see Appendix A 
The energy E is the difference between this approximation and the sum (II ip and is 
thus expected to be negligible in the limit L ^> c/T. For large and vanishing D, the 



Sudden change of the thermal contact between two quantum systems 5 





A 

4 






■ '- 








LU 




7 /'' D = ~ 




-4 


" / D = 0.02c/T 






D = 0.1c/T 




-8 


:\ D»c/T - 

■ ■ . i .... i .... i . . , . i . . . . ■ . . . . 





0.1 0.2 0.3 0.4 

TL/c 



Figure 1. Energy E as a function of length L for A = O.Olc/T and different 
values of D. E reaches a maximum at D ~ 0.03c/T. The corresponding curve is 
indistinguishable from that obtained for D = 0.02c/T. The results for D > lOOc/T 
and D < 10 _4 c/T are well-described by the large and vanishing D approximations 
(fT3f and (|T4"|) . respectively. 



wavenumbers q are equally spaced and it is then possible to evaluate E 1 see | Appendix A 
We find 

E (oo, TL/c) = | - —T 2 L sinh(2pvrTL/c)- 2 (13) 

E(0,TL/c) = - — T 2 L Vsinh(4p7rTL/c)- 2 . (14) 
c 

p>0 

Note that the energy {H)t for large D is not exactly twice the energy of an isolated cavity 
of length L, the coupling between the two cavities contributes an energy T/2. In the limit 
L ^> c/T, the energy ffl3|) reaches T/2 and (|14p vanishes. For an arbitrary characteristic 
coupling length D, E/T converges also to a finite value in this limit, as shown by Fig. (TJ 
and hence the energy of the complete system (H)t — Lc/ttA 2 + ttLT 2 /3c is essentially 
independent of D. We remark that, since the solutions of ([6]) are q ~ pit /L + 2/npD 
where p is a positive integer, for large D, E first increases with decreasing D. The 
expressions (fl3"j) and (fT4"|) are valid for any LT/c. We see that, for L c/T, the 
contribution of the coupling (jHJ) cannot be neglected as E (oo, TL/c) and Tv (0, TL/c) 
diverge differently in this limit. We obtain this diverging behavior of the energy L for 
other values of D as well, see Fig. [TJ 
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Figure 2. Equilibrium correlation function (H(x)H(x'))t as a function of x 1 — x for 
T = 5c/L, x = O.IL and 0.2L, and D = 0.1L, L and 10L. The dotted line corresponds 
to the thermodynamic limit expression (|18p . For the above parameters, the agreement 
with this approximation is excellent for x' > x. 



2.3. Equilibrium correlation functions 



It is also instructive to consider the thermal two-point correlation functions of the 
system, for example 

-A| g | 



aw)> T = £ E d£iri ^"" ] + e * w> ) < 15 > 



q=pw I L 



tanh(cg/2T) 



H — sgnlxx ) > — — — — 

8 6 v ; ^tanh(cg/2T) 



e iq(\x\-\x'\) _|_ 



Dq - 2i 
Dq + 2i 



M{\x\ + \x'\) 



where p runs over all integers and the second sum over the negative and positive solutions 
of (EJ). For vanishing and large D, this expression can be simplified by using the relation 
( lA.lj) and the function (1A.3I) . We find, for D = and D — > oo, respectively, 

(n(x)n(x')>T = - TT~ sinh ~ 2 ~{ x + ex> + ( 4 P + 1 + e ) L 



(n(x)n(x'))^ 



2c 

TtT 2 

2c 



p.c 



Q(xx') sinh -2 — yx + eo/ + 2pL 



(16) 
(17) 



where p runs over Z and e over {—1, 1}. For L ^> c/T, the only terms which contribute 
significantly to the sum (TIBl) are (e,p) = (—1,0), (1,-1) and (1,0). The last two are 
important only for (x,x r ) in the vicinity of (—L,—L) and (L,L). The vanishing of 
(I17p for xx' < shows that there is no correlation between the two cavities in the 
large D limit. The comparison of (|T6l) and ([17]) reveals that, in the large D case, the 
thermal correlation function (|T5l) is identical to that of two isolated cavities of length 
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L. Consequently, in these two limiting cases, 

7rT 2 

(U(x)U(x')) t ^ sinh- 2 [7rT(x - x')/c\ (18) 

for x and x' not too close to the cavities ends. More precisely, (U(x)U(x'))t deviates 
notably from this approximate expression if x or x' is at a distance smaller than c/T 
from —L, or L. For other values of D, the field-field correlations (|T5l) are also well 
described, in the large L limit, by (|T8l) . This can be seen as follows. For x ~ x' not 
too close to the ends — L, and L, the terms exp[iq(x + x')] and exp[zg(|a;| + \x'\)] in 
( IT5l) vary a lot from one value of q to the next and hence give vanishing contributions. 
Moreover, the normalisation factor A q is close to L~ 1//2 except possibly for the first few 
roots of ([6]), and this equation has only one solution in a given interval [pn/L, (p+1)tt/L] 
where pGN. As a result, the two remaining sums in (|T5l) are accurately approximated 
by the same integral which leads to (fT8l) . The thermal correlation function ffT5l) differs 
significantly from this approximation only for x and x' close to the cavities ends as 
shown in Fig. [2j 

For the field cd x <p, we obtain the expression 

c^ { x)8 M = f L £ [e^ - )) (19) 



8 ^ tanh(cg/2T) V Dq + 2i 



where p runs over Z and the second sum over the solutions of For D = and 
D — ► oo, this correlation function is given by fTT6l) and f|T71) . respectively, with an extra 
factor — e in the summand. For other values of D, the above arguments can be used to 
simplify (1191) . Thus, the approximate expression ( |T8i) applies to the field cd x (p as well. 
For the correlations between the fields II and d x (j), the completeness of the basis {<p q } 
implies, for any D, (d x <j)(x)U(x')) = (i/2)5'(x — x'). In summary, in the thermodynamic 
limit, the thermal correlation functions are essentially independent of the characteristic 
coupling length D. 



3. Sudden change of the coupling strength 

We consider here that the system is initially at equilibrium with temperature T and 
characteristic coupling length Dq and then evolves under the Hamiltonian ([7j) with 
D Dq. We study the time-dependent two-point correlation functions and the time 
evolution of the energy contained in one cavity. Since, for these expectation values, a 
microcanonical mixed state or a pure state of macroscopically well-defined energy are, in 
the thermodynamic limit, equivalent to a canonical ensemble [2], we evaluate canonical 
averages in the following. 
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3.1. Time-evolved field operators 



To obtain the field-field correlations at any time t, we write the time-evolved field 
operators in terms of the creation operators a\ corresponding to Dq as 

1 



d x <p(x, t) 



JHt 



d x <pe 



-mt 



Q k (x,t)al + Qt(x,t)a k (20) 



fc>0 



U(x, t) = P + iJ - e~ kA/2 Pk(x, t)a\ - Q* k (x, t)a k 

" fe>0 

The coefficients Q k and Q k are given by 

Q k (x,t) = k- 1 / 2 J2(k\q)d x ^ q 



q>0 



k 



t) = k 1 ^ 2 ~S^{k\q)(p q (x) cos(cqt) + i— sin(cqt) 



cos(cqt) + i— sin(cgt) 
■ Q 



q>0 



(21) 

(22) 
(23) 



where {k\q) = j^ L dxif k °\x)f q (x), ip q and ip k are the eigenmodes corresponding to 
D and Dq, respectively If (p k is an even function of x then only the even ip q 
contribute to the sums (12"2"|) and f[2"3"j) . Moreover, in this case, since the even eigenmodes 
do not depend on the characteristic coupling length, these expressions simplify to 
Ok = — (k/L) 1 / 2 sin(kx) exp(ickt) and Q k = {k/L) 1 / 2 cos(kx) exp(ickt). For odd <p k °\ 
it is instructive to first consider the two limiting cases of vanishing and large D. 

3.1.1. Perfectly coupled cavities For D — 0, the odd eigenmodes are <p q {x) = 
L~ l l 2 sin(gx) with the wavenumbers q = (p + 1/2)tc/L where p E N. By evaluating 
(k\q) and using the relation flA.lj) . we obtain 



Q k = 2iLA k cos(kL)k- 1/2 J^(-l) p \s{t - x - p) - 8(t + x - p) 
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_ ^± k l/2 e -ikL ^_^[T-x\ e 2ikLe{r-x) _ ^_^\T+x\ e 2ikLe{T+x) ^4) 

where p runs over all integers, r = ct/2L, x = x/2L, \t\ denotes the largest integer 
smaller than r, and e(r) = r — \ t\ . Inspection of this expression shows that the field II 
is the superposition of two waves travelling with velocity c and reflected with no phase 
shift at x = ±L. For a given x, \Q k \ takes two different values as time goes on, see Fig 
H We find 



e - lkct fl 



k 



e i(7r 2fcL)p sgn ^ cos[k(\x\ - L)) in regions A, A', 



A k Vk 

= e i (— 2fcL )(p+ 1 /2) S i n (k x ) in regions B, B', . . . (25) 

where p = in A/B, 1 in A'/B', .... The space-time regions A, B, A', B', ...do 
not depend on the wavenumber k. As we will see below, this independence plays an 
essential role in the behavior of physical properties such as the field-field correlations 
or the cavities energies. For the field d x <f), the coefficient (122]) is given by a similar 
expression. 

3.1.2. Uncoupled cavities For large D, we obtain 

n k = sgn(x)^k^ 2 e- ikL e 2lkLe{T ~ s) + e 2ikL < T+s) . (26) 

Here, the field II is the superposition of two waves reflected at x = and L (-L) in the 
right (left) cavity. The expression (125]) becomes 



e -ikct Q 



k 



e 2lkLp sgn(x) cos[k(\x\ — L)\ in regions A, A', . . . 



A k Vk 

= e~ 2ikL(p+1/2) sgn(x) cos(kx) in regions B, B', . . . (27) 

where p = in A/B, 1 in A'/B', .... Contrary to the above case, the field II is always 
discontinuous at x = 0. The field d x (j) satisfies d x <fi(0, t) = 0, i.e., the reflection is perfect 
at x = + and 0~. 

3.1.3. General case For a finite value of D, the fields II and d x <p are also superpositions 
of propagating waves but which are, contrary to the above studied special cases, both 
transmitted and reflected at x = 0. We write the coefficients flk and Ok as 



Q k (x, t) = sgn(x) (-^- - lj A k cos(kL)k^/ 2 [T+(x, t) + T^(x, t 
®k(x, t)=i(^--l)A k cos(AlL)AT 1/2 [Tt(x, t) - T k (x, t)] (29) 



EaA 2 P *<?( ci ±N) 
(30) 
k-q Dq±2i K ' 

q 

where the sum runs over the positive and negative solutions of (E])- To evaluate the 
functions T^, we first note that, since q is a solution of the summand in (!30j) can be 
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multiplied by a factor \{Dq + 2i) exp(— 2iqL)/(Dq — 2i)] m where m is any integer. We 
then rewrite Tf as 

If X(z) m+1 ze ^t±\x\)/D 
r^(x,t) = — (b dz 



nD Jc 1 - X(z) (Dk -z)(z± 2%) 



X \Dk) m+1 k e ik i ct ±\ x \) 
+ 2 \-X(Dk) Dk±2z (31) 
where X(z) = (z + 2i) exp(—2izL/D)/(z — 2i) and C is a contour enclosing the real axis. 
Expressing the above integral in terms of the residues of the integrand poles on the 
real axis leads to (ISO]) . For times t between and where t^j = =f|x|/c + 2mL/c, 
the upper (lower) part of the contour C can be closed in the upper (lower) half plane. 
Moreover, for positive m, the integrand in f[3"Tj) has poles only on the real axis and at 
z = 2%. Thus, for t > 0, the first term of the expression (!3T|) is determined by the residue 
of this last pole if m is chosen as the largest integer smaller than (ct± \x\)/2L. It can be 
shown that the functional dependence of this term on r = c(t—t^)/D is exp(— 2t)P^(t) 
where {Pm) ls a polynomial of order m — 1 (m), see Appendix B| The functions 



and 0fc change discontinuously at £ = t^. For D <^ L, they then reach the values given 
by the second term of (!3T|) in a time of order D/c. For D = 0, one retrieves the Dirac 
delta functions of (|24l) and the second term of (J3T1) gives the expression (|25l) . In the 
opposite limit, D ^> L, the first term of ( 13~TT) is negligible with respect to the second one 
which leads to the expression fl27j) . Contrary to the limiting cases of vanishing and large 
D studied above, since the successive reflections and transmissions at x = deform the 
propagating waves, and |0^| are not strictly periodic here. 

3.2. Correlation functions 

With the notations introduced in (120]) and (ED), the two-point correlation functions of 
the cavities read 

(n(M)n(^ t )> = ^ + £Reg g^g^) e^ (32) 

Since the functions 0fc and Jl^ present discontinuities at positions independent of k, as 
discussed above, the field-field correlations at a given time change abruptly at x = x m 
and x' = x m where x m = 2mL±ct, m G Z, see Fig.H]and[5j The results shown in figures 
H] and [5] are obtained by numerically solving (jSJ) and evaluating the exact expressions 
(J22J)-([23D and (j52M31l). For small D, the first term of (Ell) contributes significantly to 
the correlations (!32l) - (!34l) only if x or x' is close to x m , as shown in Fig HI This term 
becomes negligible in the large D limit, see Fig. [5j Thus, the field-field correlations are 
essentially described by the second term of (13~TT) for small and large D. Keeping only 
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Figure 4. Correlation functions (IT(x, t)H(x' , t)) (full line) and (cd x (j)(x, t)H(x' , t)) 
(dashed line) as functions of x' - x for x = 0.5L, T = 10c/ L, D = lOOi, D = O.OIX, 
and t = 0.55L/c. The discontinuities correspond to x' = ct, see Fig. [3] The spatial 
correlations of the field cd x <f) are indistinguishable from that of the field II (full line). 
The dotted line is the thermal equilibrium correlation function (| 18|) . 



this term and using the properties X(Dk)* = X(Dk) 1 = X(—Dk), we obtain 

where the second sum runs over the positive and negative solutions of (jHJ) with the 
characteristic coupling length D . The coefficients and are given by 

(a k , (3 k ) = (1, 1) for (x, a/) G A x A, A' X A',... 

= ((1 + X)/2, (1 + X*)/2) for (x, I'jGAxS.SxAA'xB,... 
= (1,X*) for (x,x') eBxB,B' xB', ... (36) 

where X = X(Dk) and the regions A, B, A', . . . are presented in Fig. [3j Therefore, the 
approximate correlation function ( 1351) . for given x and x', is periodic with period 2L/c. 
When x and x' are both in a region of type A, the expression (|35j) is identical to the 
initial condition ( fl5l) . As discussed after equation (fl8l) . the second and fourth terms of 
(j35p are negligible except for x and x' near a cavity end, and hence (1351) is practically 
equal to (Tl5l) for x and x' both in a region of type B, see Fig. HI When x and x' are 
in regions of different types, (l35i) can differ considerably from ( fl5l) . For example, for 
cavities perfectly coupled (D = 0) but initially uncoupled (D — * oo), ot k = j3 k = 
in this case and thus the spatial correlations of the field IT are reduced by a factor of 
two with respect to the initial condition ffl8]) . It is interesting to note that even for 
uncoupled cavities (D — > oo), the sudden change of the characteristic coupling length 
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Figure 5. Two-point correlation function of the field II for D = 5L, 10L and 30L. 
The other parameters are as in FigSl The curves for D > 30L are indistinguishable 
from one another. 



causes a time-evolution of the correlation functions as X = exp(—2iLk) is equal to 1 
only for D — > oo. 

For the other correlation functions, we find 

1 hp— MM / \ 

OA( M ( X >, t )) , _^ i - Ew ^ (e<*-> -e*<«'>) (37) 

+ 8c^tanh(cfc/2T) r fc6 ^ J 

where and /3& are given by (|35J) . and, for x ^ x' , 

where the sum runs over the solutions of (jHJ), the upper sign is for |x| > and 
( Tfcj 4) = (0, 0) for (x, j')eixA,5x5,i'xi',.- 

= (1 - X, X* - 1) for (x, j'jeixBjBx^i'xB',.- 
= (X - 1, 1 - X*) for (x, x') eA'xB,Bx A', A" x B f , . . . (39) 

Contrary to the correlations (135|) and (l3Tj) of a field with itself, the correlations between 
the fields II and d x (f) given by ( 1381) vanish when a; and are in the same region, see Fig. HI 
For x and x' not too close to a cavity end, (d x <p(x, t)d x i<j)(x', t)) ~ (II(x, £)II(V, t))/c 2 and 
(9 x 0(x,t)n(x',t)) ~ ±((n(x,t)n(a; , ,t)) - (U(x,0)U(x',0)))/c where the sign depends 
on the relative positions of 0, x and x'. We have seen in the previous section that there 
is no correlation between the two fields and that the correlations of a field with itself are 
essentially the same everywhere at equilibrium. Switching the characteristic coupling 



16 ^ tanh(cfc/2T) 

^ 7fce ik||*|-M| ± Ske ik(\x\+\x'\-2L) 
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Figure 6. Coupling energy as a function of time for Dq — L,T — 10c/ L and D = 0.1L 
and 0.2L. In the left part of the figure, t < to for D = 0.1L and t ~ to for D — 0.2L. 
The dotted lines correspond to the small D approximation (|43|) . 



length to another value induces correlations between the two fields and affects the whole 
system even far away from the interface between the two cavities. 

3. 3. Energy of the cavities 

We study here the time evolution of the energy content of the cavities. For that purpose, 
we first note that the energy density 

is even with respect to x. This symmetry property leads to 

(H Mt (t)) = (F right (t)) = (H)/2 - (iJ coup i ing (t))/2. (41) 

The time dependence of the cavities energies is thus simply related to that of the coupling 
energy. As the eigenmodes of the total Hamiltonian H satisfy d x tp q (CT) = d x ip q (0 + ) = 
[ip q (0 + ) — ip q (CT)] /D, the time-evolved coupling Hamiltonian Qj is proportional to 
d x <p(0, t) and hence can be expressed in terms of the functions ( |22|) . We obtain 

(iWng W> -Z D )L tanh( C fc/2r) • (42) 

The even eigenmodes do not contribute to this sum since their derivatives vanish at 
x = 0. For the odd eigenmodes, we have seen above that, for D L, the first term of 
the expression (|3T|) contributes only for t close to 2mL/c where m is an integer. The 
second term of ( 1311) gives, using \X(Dk)\ = 1, 

A 2 ke~ kA 1 

<iW., s «> t ^(ck/2T)4 + (Dkr < 43 » 
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Figure 7. Coefficient O qq / as a function of frequency u> qq > = c(q—q') for the correlation 
function (II(x, t)U(x',t)) with x = 0.5L, x' = 0.51L, D Q = 0, D = 0.1L and T = 10c/ L. 



where the sum runs over the wavenumbers k corresponding to odd <p^ . For a small 
but finite D, the coupling energy deviates from this value only for short periods of time 
when t is much smaller than a characteristic time to which increases with decreasing D. 
In this short time regime, the energy density (flOl) is essentially constant except in two 
regions wich propagate with velocity c and are reflected and transmitted at x = 0, L and 
—L. The multiple peak structure shown in left part of Fig. [6] results from the preceding 
reflections and transmissions at x — 0. The temporal evolution of (-f/coupimg(0) becomes 
more complicated for t > to, see Fig. [6j For D ^> L, the contributions of the first 
and second terms of fl3T|) to the functions 0^(0, t) are of the same order and hence 
the approximate expression (J15|) does not apply. For D — > oo, d x (f>(0, t) vanishes, as 
mentioned after ( 1271) . However, because of the factor D in ( H2|) . some care must be 
taken in the evaluation of (H coupling (t)) in this limit. For large D, the wavenumbers q 
are practically the multiples of tc/L except the lowest one q ~ (2/DL) 1 ' 2 . The terms 
of the sum (122]) are thus essentially equal to their infinite D values except the first one 
and 

itt tj.w ■ i( ct \ c ^ A 2 k sin 2 (kL)e~ kA 

(H^t)) = s,,^ (^-^ j - g ,\ anh ' (c ,; 2r) (44) 

varies sinusoidally with a finite amplitude. 



4. Comparison to perfect gas 

The time- dependent average energies (Hi e { t {t)), (if rig ht(t)) and (H conp \i ng (t)) , and the 
correlation functions (!32~ l) -( !34l) can be written as 

(6(t)} = J2°<M' eit{eq ~ eq,) = J due itu> 0{u) (45) 

9,9' 
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where O stands for H\ eit , . . . , the sums run over the positive and negative q corresponding 
to the characteristic coupling length D, O qq i are appropriate coefficients, e q = cq and 
= J2 qg >>o O qq '5(uj — e q + e q /). We remark that the even eigenmodes <p q contribute 
only to the steady component of (0(t)), see the discussion after equation (|23l) . The 
expectation values of the single-particle observables of an isolated perfect gas, its particle 
number density at a given point for instance, are also given by expressions of the form 
( 145]) with the single-particle energies e q [10] . Consider as an example a one-dimensional 
perfect gas confined in a box of length L. In this case, the single-particle energies are 
e q = p 2 7r 2 /2mL 2 where m is the particle mass and p an integer. The frequencies e q — e q > 
are thus regularly spaced and (0(t)) is periodic with period 4mL 2 /ii. However, in the 
Joule expansion studied in Ref. [10] , for times t <C mL 2 , the function O in the integral 
expression f]45]) of the gas density profile, can be approximated as a continuous function 
plus a term O^diu) and (0(t)) essentially relaxes from its initial value to Oqo. 

The situation is radically different for the free field system considered in this paper. 
There obviously also exists a 5(cj) contribution to 0(u) but {Of)) does not relax. First 
of all, we observe that, in the limiting cases of vanishing and large D, the function O is 
a sum of equally spaced Dirac delta functions and {Of)) is strictly periodic with period 
2L/c as shown by ( 125]) and ( 127]) . For a finite D, the frequencies u qq i = t q — t q i = c(q — q') 
are not exactly equal to multiples of nc/L but most of them are close to these values, see 
Fig. [7] As it is clear from equation ([6]), the spacing between consecutive frequencies u qq > 
varies and goes to zero for u qq i — > pnc/L where p G Z. We see on Fig. [TJthat O qq i — > 
in this limit. Time regimes may then exist in which 0(uj) can be approximated as a 
sum of smooth functions with finite support. However, these functions clearly vanish 
at the right end of their support interval but not at the left one. Consequently, their 
Fourier transforms do not vanish fast enough at long times to lead to a relaxation 
behavior of {0(f)). We find similar functions O for the energy {H couv \ irig f)) with the 
noticeable difference that the coefficients O qq i corresponding to the lowest frequencies 
uj qq i ~ ±2c(2/ DL) l l 2 are dominant in the large D limit, giving rise to the sinusoidal 
behavior discussed at the end of the previous section. 

5. Conclusion 

In this paper, we have studied how changing the thermal contact between two quantum 
systems whose elementary excitations are noninteracting bosons, affects these systems. 
More precisely, our model consists of two one-dimensional cavities of equal length 
coupled to each other by one of their ends. The coupling we have considered is 
thermodynamically small, as shown by the fact that the equilibrium total energy and 
correlation functions are practically independent of its strength, and allows to describe 
from perfectly coupled to completely uncoupled systems. We have seen that a sudden 
change of this thermal contact has a considerable impact on the correlation functions 
of the two cavities. They do not return to equilibrium but essentially oscillate with 
a period equal to twice the time required for a signal to propagate from one end of a 
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cavity to the other. Their exact time evolution is mainly controlled by the final coupling 
strength which determines the reflection and transmission coefficients at the interface 
between the two systems. We found a similar behavior for the energy content of the 
cavities. For weak coupling between them, it varies sinusoidally with a period set by 
the coupling strength. The absence of relaxation behavior is not simply a consequence 
of the fact that the considered Hamiltonian is quadratic. The Joule expansion of a 
perfect quantum gas has, for example, been obtained in Ref. [TOJ. As discussed in 
the previous section, this difference in behavior comes from qualitative differences in 
the elementary excitation energy spectrum. Clearly, the simple dispersion relation 
of our one-dimensional continuous field model, plays a crucial role in the temporal 
quasiperiodicity obtained. It would be thus interesting to study higher- dimensional and 
lattice free field systems. Are interactions between the bosonic excitations necessary to 
ensure the stability of the thermal equilibrium ? 



Appendix A. Coupling energy in limiting cases 

To evaluate the energy (fTTjl for vanishing and large D, we make use of 

J2F(p)e ipx = J2 J dqF(q)e iq{x - 2pn) (A.l) 
p p 
where p runs over all integers and F is any function. 
With this relation, we obtain, for large D, 

(H) T = - + —LT 2 G(2TA/c, 0) + —LT 2 V G(2TA/c, ApTL/c) (A.2) 

2 TTC 7TC Z — ' 

p>0 

where 

f°° k 

G(e,x) = Re dk — e - ke+ikx . 

J q tann k 

Using the expansion tanh -1 k = 1 + 2 J2 p>0 exp(— 2kp) and ( 1A.1I) . we find, in the limit 
e -> 0, G{e,x) = -7r 2 [2sinh(7rx/2)]- 2 for x ^ 0, and G{e, 0) = e~ 2 + tt 2 /12. Replacing 
these expressions into (1A.2|) leads to (fT2]) with (TT3|) . 

For D = 0, the wavenumbers of the odd eigenmodes are q = (p + l/2)ir/L whereas 
q = pir/L for the even ones. The contributions of the odd and even eigenmodes are then 
different. We obtain 

E(0, cT/L) = -—LT 2 Yu p - —LT 2 V (-l)X (A.4) 

7TC TTC 

p>0 p>0 

where u p = sinh(2p7rTL/c)~ 2 , which gives ([T 



Appendix B. Time-dependence of the field operators 

The polynomials and P~ defined after (jHTl) are given by 

Pl(r) = P m -i(r) = ^ B j c dzF{z) (B.l) 
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= ze«* ( z + 2f\ m j. _ 2i2i/I3 _ _ , -i 
v ; Dk- z \z-2ij lK ' J 

The above integral is determined by the residue of the function F at z = 2i. We write 

m + e) = V + + [£e (-4, i/D _ 4i _ e] - 

^ p! DA; - 2« - e L J 

p>0 ^ 

(ir) p 



p,r>0 1 



-A r e p+r - m . (B.2) 



With these notations, fIB.ll) becomes 

p=0 

for m > 1 (P + = 0). We find, for example, 

P+(t) = - —— (B.4) 

D Dk — 2i K J 



16 1 



iD Dk - 2i 



1 i 



-■ ~ 7 ( e " 4L/D + 3 ) + 



Dk-2i 4 



(B.5) 
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